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SUMMARY 


During the first few seconds of the space shuttle trajectory, the solid 
rocket boosters will be in the proximity of the launch pad. Because of the 
launch pad structures and the surface of the Earth, the turbulent mixing 
experienced by the exhaust gases will be greatly increased over that for the 
free- flight situation. In addition, a system will be present, designed to pro- 
tect the lifting vehicle from launch structure vibrations, which will inject 
large quantities of liquid water into the hot plume. This report studies the 
effects of these two phenomena on the temperatures, chemical composition, and 
flow field present in the afterburning solid rocket motor exhaust plumes of the 
space shuttle. The study includes results from both a computational model of 
the afterburning and supporting measurements from Titan III exhaust plumes 
taken at Kennedy Space Center with infrared scanning radiometers. 


INTRODUCTION 

As part of a continuing assessment of the environmental effects of space 
shuttle operations, the National Aeronautics and Space Administration has been 
studying the formation and behavior of rocket exhaust clouds from solid pro- 
pellant rocket motors. Although information is readily available as to the 
composition of the rocket exhaust as it leaves the motor’s nozzle exit plane 
(see, for example, ref. 1), significant changes occur in the composition of the 
plume as these hot exhaust products afterburn in the atmosphere. An understand- 
ing of these afterburning processes is necessary so that a reliable source term 
can be formed which will serve as input to diffusion models and form a basis 
for determining which products are of ecological significance. The purpose of 
this study is to predict a source term for those exhaust products which are 
released in the proximity of the launch pad. 

From free-flight calculations it is known that the solid propellant motors 
of the shuttle will release substantial quantities of such gases as HCl and NO 
into the troposphere and that material released below the local inversion 
layer will rise to form a stabilized exhaust cloud. Approximately two-thirds 
of this material is released in plumes whose afterburning properties can be said 
to occur under free-flight conditions. The remaining material is released close 
enough to the ground and launch pad structures to necessitate a separate con- 
sideration of afterburning with boundary conditions distinct from those of free 
flight. 

Studies of the afterburning processes occurring within the plumes of the 
Titan III and shuttle solid rocket motors in free flight were reported in ref- 
erences 2 and 3. These studies included temperature and velocity distributions 
within the plume as well as an analysis of the major chemical constituents. 

The altitude range considered in these studies was from 0.7 to 15 km. Below 
an altitude of roughly 0.5 km, free-flight turbulent mixing rates and ordi- 
nary atmospheric boundary conditions can no longer be considered valid. The 



interaction of the exhaust plume with the launch pad structures, the flame 
trench and associated "deluge" water, and the Earth's surface may all have 
major effects on afterburning. 

The differences between afterburning near the launch pad and in free 
flight that are considered in this study can be classified into two main cate- 
gories. The first is that the rate of turbulent mixing of air into the plume 
is significantly increased as a result of the gases which flow out of the noz- 
zle and hit the pad deflectors and structures. The other significant occur- 
rence which affects only early launch afterburning of the space shuttle is the 
initiation at a few seconds before lift-off of the so-called "acoustic water 
suppression system." This acoustic water system is described in figure 1. The 
basic purpose of the system is to protect the vehicle from vibrations in the 
launch pad structures caused by strong vibrations set up by the motors of the 
shuttle boosters. These vibrations are to be suppressed by injecting large 
quantities of water into the gas flow from the motor. This injection will 
occur near the firing pad and in the flame trench. The flame trench is a large 
open pit reinforced by concrete which runs under the firing pad and which chan- 
nels the afterburning gases away from the launch area. For the space shuttle, 
the solid rocket motor gases will travel a path through the trench which is 
about 80 m in length. 

Present plans are that at 4 seconds before launch water will be injected 
under the firing pad and into the flame trench at a rate of 16 000 il/sec on the 
solid motor side (fig. 1) until at least 8 seconds after lift-off. This rate 
represents a mass flow rate of about 16 000 kg/sec from lift-off through this 
period; and a large percentage of this flow will be entrained into the plume. 

The mass flow rate of the plume is somewhat less than this rate, being roughly 
10 000 kg/sec. 

Although the effect of water injection on afterburning occurs primarily 
through temperature quenching, the effect is also closely coupled to the tur- 
bulent mixing processes. Since modifications to the basic physical model are 
required to account for this two-phase flow, the water injection problem was 
studied first. The following sections describe the basic numerical plume model 
and the development of a simple kinetic rate model for water injection. Results 
showing the effects of entrained water on afterburning are presented along with 
an assessment of the accuracy of the water injection model. Correlations are 
then made with scanning infrared radiometer photographs of actual Titan III 
plumes to determine the magnitude and effects of the enhanced turbulent mixing 
based on a simple mixing length model. 


SYMBOLS 


A area of plume 

a speed of sound 

C rate of radius of droplet reduction 

D diffusion coefficient of water in air near droplet 
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function 

thermal conductivity 

forward reaction rate 

heat of vaporization of water 

arbitrary third body in a three body reduction, also Mach number 
mass 

Avogadro's number 

pressure at droplet surface 

pressure far from droplet surface 

radius of nozzle 

universal gas constant 

radius of water droplet 

radius of droplet at final time tf 

radius of droplet at time t = 0 

temperature 

time 

final time 

velocity of plume gases 

molecular weight of mixture of exhaust gases and air 
mole fraction; subscript denotes compound 
distance measured radially 
scaling factor for viscosity 
cofficient of viscosity 

density of mixture of exhaust gases and air 
density of water at droplet surface 
density of water far from droplet surface 
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Subscripts: 


c inner core 

dv developed region of mixing 

r at radius of droplet 

oo far beyond radius of droplet 


MODEL DEVELOPMENT 
Basic Numerical Model 

The basic numerical model used to simulate the afterburning processes is 
the low-altitude plume program (LAPP) code. (See ref. 4.) This is a steady- 
state axisymmetric code which couples the effect of kinetic rate chemistry to 
turbulently mixing fluid flows. The code uses an explicit-implicit finite- 
differencing scheme to "march” downstream of the nozzle exit plane. 

As inputs to the code, it is necessary to specify the conditions at the 
nozzle exit plane as well as the conditions of the unperturbed atmosphere into 
which the nozzle gases are flowing. The number of species in the plume, the 
chemical reactions which they undergo, and the reaction rate coefficients must 
be included. It is also necessary to specify the rate of turbulent mixing. 

The LAPP code contains several options for turbulent mixing models. The 
model used in this study is the Donaldson-Gray model (ref. 5) which was chosen 
because of its successful use in prior studies (ref. 6). The Donaldson-Gray 
model contains an undetermined mixing coefficient which is usually set equal 
to unity under tropospheric free- flight conditions. This parameter is referred 
to as ot. Later in this study, a reasonable value for this coefficient in the 
highly turbulent flow near the launch pad was determined. 

As discussed in the "Introduction," water is injected into structures 
beneath the pad at a rate of 16 000 kg/sec starting at 4 seconds prior to 
launch. When the actual firing of the solid motors occurs, these gases will 
essentially be flowing into a "lake" of water. Much of the water will be 
entrained into the plume which, in turn, will vaporize the water into steam. 
This highly endothermic vaporization reaction (2261 J/g) will tend to quench 
some of the high temperature afterburning. 

In figure 2 the importance of the phase change is shown for the prelimi- 
nary case of a plume whose exit plane constituents are precisely those of a 
shuttle SRM, but where no chemical reactions are allowed to take place. In 
this inert plume, all cooling takes place by the mixing of a cool substance 
with a warm substance and only nonchemical heat transfer takes place. In the 
calculations leading up to this figure, the phase change was included by using 
an early thermodynamic equilibrium model. In this figure is plotted the tem- 
perature along the plume center line as a function of distance downstream of 
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the nozzle exit plane. As can be seen, where the phase change was present 
(with liquid water in the atmosphere), the cooling took place significantly 
faster than with air only or air plus water vapor. It was thus considered 
important to treat the phase change in any set of detailed calculations. 

The LAPP code makes the assumption that all substances in the plume are 
ideal gases. Clearly, for the study under consideration where liquid water is 
vaporized, such a single phase model is inadequate. After investigating several 
alternate approaches, it was decided to include the effects of the acoustic sup- 
pression water by having a mixture of air and liquid water entrained into the 
plume in the first 80 m (the path length of the flow field from the nozzle at 
the pad through the flame trench) and entraining pure air from 80 m until the 
problem terminates about 1 km downstream of the exit plane. 

To account for the phase change during evaporation, liquid water is con- 
sidered a new chemical species which obeys the equation of state of a perfect 
gas as required by the code, but with thermodynamic properties (heat capacity, 
specific enthalpy, etc.) appropriate to liquid water below 373 K and appro- 
priate to water vapor above 373 K. The heat of vaporization is included by 
means of an imaginary chemical reaction: 

Liquid water + M — ► H 2 O + M (la) 

irreversibly 

This reaction is endothermic, 40.70 kJ/mole, which is the heat of vapor- 
ization of a mole of liquid water. M is an arbitrary third body. 

It was found during the course of the investigation that because of the 
preponderance and stability of the N 2 molecule. 

Liquid water + N 2 — ► H 2 O + N 2 (1b) 

irreversibly 

could be used with negligible differences in the results. 


Rate Coefficient Calculation 

In addition to specifying the constituents of a reaction, it is also 
necessary to denote the value of the reaction rate coefficient appropriate to 
the reaction. The reaction rate coefficient appropriate to equation (lb) is 
calculated as follows: 


Let Xiiquid water equal the mole fraction of liquid water and Xjjg equal the 

mole fraction of N 2 . From the formulation of the chemistry in the LAPP code 
(ref. 4), it can be deduced that 


^liquid water^^) = ^liquid water^T = 0) exp 


Pmix 

-NAkf(T)— XNgtl 
"mix 
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where is Avogadro’s number, kf(T) is the unknown forward reaction rate 

coefficient dependent upon temperature T, t is time, Pmix density 

of the mixture, and molecular weight of the mixture of plume 

gases, air, and water. 

Throughout the wide range of physical and thermodynamic conditions encoun- 
tered in this study, it was invariably found that the overwhelming preponder- 
ence of liquid to vapor transition occurred in the stream tube most distant from 
the center line with a temperature above 373 K. It was further observed that 
many quantities in the stream tube were relatively constant. For example, 

Wmix stays within 10 percent of 27 g/mole, Pmix stays within 10 percent of 
1 mg/cm3, and stays within 20 percent of 0.8. These results hold for all 

values of downstream distance. Using these values simplifies the rate expres- 
sion to 

^liquid water^^^ ~ ^liquid water^"^ ~ exp{-kf>(T) 1.8 x 10^^ t} 


Since more than 98 percent of all vaporization takes place in the outer- 
most stream tube at a temperature of approximately 373 K, it is not surprising 
that once the value of kf is fixed in this stream tube, the behavior of the 
plume is relatively insensitive to the temperature dependence of kf in the 
higher temperature stream tubes. Several different forms of kf(T) were tried 
and included: 

kf(T) = kf(373 K)^ (2) 


kf(T) = kf(373 K)T (3) 

kf(T) = kf(373 K)t2 (H) 

Using equations (3) and (4) was found to be somewhat more expensive compu- 
tationally than using equation (2) and did not affect plume behavior signifi- 
cantly. Hence, kf(T) = kf(373 K)\|T (eq. (2)) was used throughout the study. 


By using the form of the rate coefficient described and assuming a value 
for kf(373 K) , the model predicts a specific ”e-folding" time (or time to 
reach 1/e times its original value) for water evaporation at every point in the 
flow field. Examples of these times are given in table I. To determine the 
proper value of kf(373 K) , it is necessary to have some idea of what the true 
e-folding time is for water droplets evaporating in this situation. This ques- 
tion is treated in the next section. 


Heat and Mass Transfer in Water Droplets 

In this section an attempt is made to predict the rate of evaporation of 
water droplets. This study follows closely the techniques developed in refer- 
ence 7. The basic assumptions in this calculation are that (1) the heat flow- 
ing into the droplet from the carrier gas is equal to the heat lost by the 
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droplet through evaporation, and (2) the mass flow from the droplet can be 
viewed as a diffusion process. 

The heat transferred into a droplet of radius r with thermal conduc- 
tivity K, and a surface temperature Tp from a carrier gas with temperature 
Tgg will be given by 4iTKr(T^ - T^). Equating this relation to the heat lost 
by the droplet through evaporation 


-L — = 4irKr(T - T^) (6) 

^4- OO i 


where L is the heat of vaporization of the water. 

If water molecules are thought of as leaving the surface as a result of a 
diffusion process, then 

dm 

= 4TTrD(pr, - pj (7) 

dt r oo 


where D is the diffusion coefficient of water in air near the droplet, Pj, 
is the density of water at the surface of the droplet, and p^^^ is the density 
of water at a distance large compared with r. 

This model neglects the phenomenon of surface tension and possible spall- 
ation of the droplet. Since the LAPP code assumes all particles and gases flow 
at the same speed in each stream tube, conduction effects are also neglected. 
The technique just described is most appropriate for a quiescent moderate tem- 
perature medium rather than the hot turbulent edge of a rocket plume. Never- 
theless, since diffusion and heat transfer within the droplet are dominant 
processes, the resultant e-folding time should be reasonable. 

Equating dm/dt in equations (6) and (7) and using the ideal gas equation 
to define Pp and p^ results in the following equation: 


K(T„ - T^,) 

£ = D(Pr - Pj 




( 8 ) 


where 


K = 2.6 X 10“^ J/cm-K-sec 


Rg = Gas constant 
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L = 2.26 X 103 j/g 


D = 0.219 cm^/sec 


The quantity P„/T^ is assumed small compared with Pp/Tp and the term is 
dropped. 

The undetermined quantity is Tp or more appropriately Tp - T^^, the 
driving force behind the heat transfer. 

By assuming that Pp is the saturation vapor pressure of water at the 
droplet temperature Tp, and by consulting appropriate vapor pressure tables, 
T~ - T^ is found to be about S5 K. 

i oo 

By assuming a particle density of 1 g/cm3 for water, equation (6) becomes 


-Lr — = k(T„ - Tp) (9) 

dt 

Under the assumptions that L, K, and T^ - Tp are constant with respect to 
time, integration of equation (9) implies that 

r2 = Pq 2 _ ct (10) 

where 

2k(T„ - Tp) 

C = (11) 

L 


In this case 


C = 1956 (jim)2/sec 


As a working definition of evaporation, a droplet will be said to have 
evaporated when its volume has reached approximately 1 percent of its original 
volume. Since e“^ = 0.01, this condition corresponds to 4 e-folding times. 

In terms of radii 

Tq = rf = 4.65rf (12) 


Equations (10) to (12) lead to the conclusion that the rate of evapora- 
tion depends on the initial size of the droplet. In table II is presented the 
evaporation time as calculated from these equations as a function of initial 
particle radius. 
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Water droplets from spray nozzles such as the ones used in the acoustic 
suppression system generally range anywhere from 10 pm to 0.5 cm. As shown 
in table II, this value corresponds to evaporation times from 50 msec to 
extremely long times. In the system at Kennedy Space Center the large end of 
the droplet size spectrum will probably be dominant. By comparing the e-folding 
times in table I with the values in table II, the appropriate rate coefficients 
are deduced as 10"^®\)f^ or lower. 


These results must be used with a great deal of caution. For particles 
larger than a few hundred micrometers, the moving of mass away from the water 
droplet by diffusion is a severe limitation. The medium in which the particles 
move is in turbulent motion, large thermal gradients exist in the regions under 
consideration - far larger than the difference between droplet temperatures and 
373 K - and collisions occur. All these effects are neglected in the heat flow 
and diffusion calculation made. 


For this reason, error bands of a factor of 103 are put on the value of kf«, 
and the water entrainment effects were studied with kf(373 K) as high as 10”^^ 
and as low as 10“20_ Most attempts to run the calculation with kf(373 K) higher 
than 10”^^ were frustrated by the prohibitive amount of computer time necessary 
to run these faster cases. Since large droplet sizes are anticipated, the lower 
rate coefficients will apply to most droplets. 


MODELING RESULTS 
Water Entrainment Calculation 

Four cases of acoustic water injection were studied. The first was an 
excess water flow condition which is appropriate for the initial firing where 
there is a large excess of water in the flame trench. Here a mass flow of 
water equal to twice the mass flow of the exhaust gases is entrained into the 
plume in the first 20 m. Another case studied was the matched flow condition 
where the amount of water entrained was equal in mass flow to the exhaust gas 
flow. This is roughly appropriate for a few seconds after lift-off. Cases for 
one-half and one-quarter of the mass flow of the gas were also studied. In every 
case an additional one-third of the value of mass flow is entrained from 20 to 
80 m because of additional water injection. 

Among the primary effects of the suppression of afterburning as a result 
of the quenching water is that increased amounts of carbon monoxide survive the 
afterburning process over the free-flight case. To obtain an estimate of this 
increase that would tend to be conservative from an environmental or safety 
standpoint, the high value of kf(373 K) = 10“^^ was used in a detailed set 
of calculations. The chemical reaction scheme besides that of the water evap- 
oration reaction consists of 22 reactions involving 18 chemical species. These 
reactions and their rate coefficients are presented in table III. 

The nozzle exit plane conditions used in this study are given in table IV. 

An entrainment coefficient of a = 1 was used in this study which, as will be 
seen, slightly underestimates the rate of entrainment of water, but does not 
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cause great errors in the final results. The results are shown in two forms. 

The form of major importance introduces the concept of number of grams of 
species per downstream meter of plume. As the gases flow downstream of the 
exit plane, they entrain cool air and water and they react chemically. As 
the entrainment process proceeds downstream, the gases rapidly cool to the 
ambient temperature and slow to the speed of the air (in the frame of refer- 
ence of the motor). As the gases cool, the kinetic reactions slow considerably 
so that several hundred meters downstream of the exit plane, the chemical 
reactions are proceeding at a rate that is negligible in comparison with the 
rate at the exit plane. If one were to take a cross-sectional slab of the 
plume in the radial direction with a thickness of 1 m in the downstream direc- 
tion at some distance downstream, such a volume would contain so many grams of 
HCl, so many grams of CO2, etc. The number of grams of such a species per 
downstream meter of plume gives one an estimate of the quantity of species being 
deposited in the cool, slowly reacting region of the rocket plume. 

Table V is such a list of species for a distance of 1 km downstream of the 
exit plane for the cases delineated. There are several things to notice in this 
table. The amount of carbon monoxide in the plume that does survive afterburn- 
ing gets larger as the amount of entrained water increases. However, for the 
mass flows of water considered, at no time are there particularly large quan- 
tities of carbon monoxide, at least when using kf (373 K) = 10 “^°. Another 
point to recognize from this table is that the amount of AI2O3 fluctuates 
slightly. In the model of afterburning being considered here, AI2O3 is an 
inert quantity. Hence, the quantity 

(pUA);^l203 = Constant 

That is, its mass flow must remain constant. If, as one would expect, at 1 km 
from the exit plane, the speed of the aluminum oxide particles is roughly the 
same for all cases of water entrainment, then the amount of AI2O3 in g/m must 
also be approximately the same. The table bears this conjecture out, the small 
relative variations in AI2O3 being due more to "noise" in the numerical scheme 
than true speed variations. (See ref. 3 .) The small percentage variations in 
the quantity of CO2 may also be considered numerical "noise" variations except 
for that effect of the nonconversion of CO to CO2 during suppressed afterburning. 

On the other hand, the large relative drop in the production of nitric 
oxide with increasing amounts of deluge water is due to the quenching of after- 
burning. Similarly, the production of CI2 from HCl falls off for the same 
reason. The success of these large amounts of water in suppressing the high 
temperature region is shown in figure 3 . This figure shows center-line temper- 
ature as a function of downstream distance for the various cases. The temper- 
ature never exceeds that of the nozzle exit plane for the excess flow case; 
whereas with unquenched afterburning, temperatures rise several hundred degrees 
above that of the exit plane. 

In figure 4 is shown the chlorine partitioning of the plume among the 
various possible species present along the center line of the plume. The 
decrease of all concentrations of chlorine-containing species at large distances 
downstream of the exit plane is due to dilution rather than to afterburning. 
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In figure 5 a similar plot is made for combustion species. Here are 
plotted the mole fractions of H2, CO, H2O, liquid water, CO2, and OH in the 
plume as a function of the distance downstream from the nozzle exit plane. As 
can be seen, as afterburning proceeds, the initially dominant concentration of 
CO drops as its molecules move downstream and are oxidized to become CO2. The 
small quantities of CO which appear at large distances are entrained from those 
regions of the plume in which CO is not completely burned. As can be seen from 
this figure, the concentrations of H and OH drop off rapidly as they combine 
to form water vapor. In figures 6 and 7 the same chemical species are shown 
in g/m. (Liquid water is not shown in this graph because its value quickly 
becomes too large.) The asymptotic values of these curves represent the quan- 
tities of species deposited in the low temperature regions of the pliune. 

It should be noted also that in the various cases considered, larger and 
larger amounts of liquid water are vaporized. Consequently, more and more of 
the energy of the plume is used to vaporize the water; as a result, the amount 
of energy available for cloud buoyancy is reduced. 


Sensitivity to Water Evaporation Rate Coefficient 

In addition to the results obtained, a parametric study was undertaken 
to determine the sensitivity of these results to uncertainties in the evapora- 
tion reaction rate coefficient. The excess water flow situation was used as 
a test case. In this study the chemical scheme used was identical with that 
of table III except that the reaction rate coefficient of reaction (20) was 
varied as shown in table VI. The results of this study are also presented in 
this table. Again, as can be seen in the variation of AI2O3, variations of 
3 percent or less must be considered physically meaningless. However, several 
overall trends can be ascertained from the results. As the rate coefficient 
for evaporation of water becomes larger, the amount of HCl increases at the 
expense of CI2 ; the amount of carbon monoxide that survives afterburning 
increases significantly; and the production of nitric oxide drops rapidly. 

The temperature quenching resulting from the various water evaporation rate 
coefficients is illustrated in figure 8 by the center-line temperature varia- 
tion along the plume. As can be seen, the larger rate coefficient reactions 
suppress afterburning strongly. 

In figure 9 the quantity of chlorine-containing species (HCl and CI2) is 
plotted as a function of downstream distance for several water evaporation 
coefficients. As can be seen, the quantity of CI2 drops dramatically as the 
rate coefficient increases whereas HCl quantities remain essentially unchanged. 
In figure 10 the quantity of NO produced is seen to drop rapidly with increas- 
ing rate coefficient whereas CO production rises (fig. 11). 


FIELD MEASUREMENTS OF PLUME TURBULENCE NEAR THE LAUNCH PAD 

Another effect which the launch pad structures will have on low-altitude 
afterburning is to greatly increase the turbulent mixing occurring in the flow 
field over the free-flight mixing. Since this increase in turbulence will have 
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directly measurable effects on plume shape and temperatures, measurements of 
these parameters have been made with infrared scanners to determine the amount 
of mixing. 

Infrared scanning radiometry of rocket exhaust clouds has become a standard 
part of the NASA launch vehicle effluents measurement techniques. Several com- 
mercially available radiometers sensitive to infrared radiation from 3 to 5.6 pm 
have been used to study the formation and diffusion of rocket exhaust clouds. 

The details of this measurement technique can be found in reference 8. 

Much like a camera sensitive to visible radiation is able to form a tele- 
vision picture of moving objects at various radiant intensities, so can the out- 
put of a camera sensitive to infrared (heat) radiation be converted to images 
on a cathode ray tube (CRT). By studying the thermal heat patterns of the plume 
of a lifting rocket vehicle, it is possible to directly infer both mass and 
temperature distributions in the plume. 

The heat pattern of a solid rocket motor plume which is in free flight 
is shown in figure 12. This infrared scan was taken of a Titan vehicle on 
March 14, 1976. The altitude of the vehicle is roughly 3 km. The horizontal 
white line in the picture is a marker which can measure intensities. The 
isotherm shown in the picture (the white vertical line in the center) is 
roughly an indication of 15 K above ambient or higher. By measuring the 
intensities along this plume, the isotherm indicating 100 K above ambient is 
present in the thick part at the top of the picture and is well within the thin 
vertical plume shown in the picture terminating about 1 km behind the nozzle 
exit plane. 

In contrast to this free-flight behavior of a constantly narrowing 100 K 
above ambient isotherm, figure 13 was taken of a lifting Titan vehicle on 
May 20, 1975. The altitude of the vehicle is about 1.1 km. The plume shown 
is entirely the 100 K above ambient or higher isotherm. As can be seen, 
although the part of the plume from 1.1 km to near the pad is constant in 
width or even decreasing, once the plume hits the ground a noticeable spread- 
ing of this high temperature region occurs with a fairly abrupt termination not 
too far downstream of the area where the spreading first occurs. 

This is shown in figure 14 in normal and downstream coordinates. Such 
behavior can be simulated on the computer by assuming more rapid turbulent 
mixing . 

With the after-deflection flow field data such as those in figures 13 and 
14 as guides, an analytical study of increased turbulence at the pad was under- 
taken. No attempt was made to develop a new theory of turbulent flow in this 
kind of situation. Several empirical models based on the concept of mixing 
length are available at the present time in LAPP. A reasonable formalism for 
the situation under study is given by Donaldson and Gray in reference 5. This 
model divides the flow field into an inner core flow and an outer developed 
flow. In the core region, the coefficient of viscosity 

p = afg(M,a,r,p ,u) 
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where fg is a function of local Mach number M, speed of sound a, radius of 
the region r, density p, and velocity u. In the developed region 

M = otf(jv(M,a,r,p,u) 

where f^jy is a different function of the same variables. Thus a is a cru- 
cial parameter since it scales the mixing rate in both regions. The param- 
eter a is set equal to unity in the free-flight case; for greater turbulence 
it attains higher values. 

For the cases studied, the vehicle is on the launch pad and the nozzle 
exit plane is in the proximity of the flame trench. The turbulent mixing param- 
eter a is assumed to be unity within the trench. From this point on until more 
than 1 km downstream of the exit plane, several different values of a were 
used in a parametric series of runs which included a = 1, 5, 10, 20, and 100. 

The values 20 and 100 are physically unrealistic in all but the most violent 
of mixing situations. 

Consider now figure 15 in which is plotted the 400 K isotherm for the vari- 
ous cases of a used. Note that the use of a = 5 from 80 m onward gives a 
fairly good representation of the measured isotherm (see figs. 13 and 14) in 
the downstream coordinate system. Based on this turbulence correlation with 
the infrared results for the full effects of the launch pad complex on after- 
burning, the most reliable results are obtained by using a = 5. 

The sensitivity of the center-line temperature to these turbulence varia- 
tions is shown in figure 16. In this figure is plotted the center-line temper- 
ature of the plume as a function of distance downstream for several values of 
a. As can be seen, the more violent the turbulent mixing, the more quickly 
afterburning is terminated. The quenching for a = 5 is seen to be much more 
rapid that for the free-flight case (a n 1 ) . 

In table VII is shown the asymptotic values of the amount in g/m of vari- 
ous species produced. Again this condition represents the excess water case 
with kf(T) = lO-l^^/f'. As a increases, CO and CI2 production increases and NO 
production seems to fall. All values have been normalized so that at a down- 
stream distance of 1 km, the value of AI2O3 in g/m is the same in all cases. 

In table VIII the combined effects of increased turbulence and excess 
water flow are contrasted with what would happen under the same nozzle exit 
plane conditions for the free-flight case. 


CONCLUDING REMARKS 

A computational model of the afterburning processes occurring in a space- 
shuttle-like solid rocket motor plume has been used to study the species com- 
position and flow properties of the plume in the proximity of the launch pad. 
The model treats parallel turbulent mixing coupled with finite-rate chemistry 
between the plume and the ambient atmosphere. The principal perturbations 
considered are the temperature quenching due to injection of large quantities 
of water into the launch pad flame trench and the increased turbulence due to 
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both the injection of water and the proximity of the plume to the launch pad 
structures and the surface of the Earth. 

The evaporation of the injected liquid water has been modeled as a kinetic 
rate process. A relation between the effective evaporation rate coefficient 
and the water droplet size has been estimated from heat and mass transport con- 
siderations. Afterburning calculations were then performed by using rate coef- 
ficients corresponding to droplet radii between 0.1 pm and slightly less than 
100 pm for amounts of water entrained varying from 0.65 to 2.4 times the exit 
plume mass flow. Radii greater than these values are also considered. In all 
cases the predicted plume temperatures are significantly lower than those in 
the free-flight case due to the evaporation of water. 

The increased turbulence near the launch pad has been modeled by varying 
the eddy viscosity coefficient in the turbulent mixing models. Correlations 
were then made between the predicted 100 K above ambient isotherm and that 
measured with a scanning infrared radiometer of a Titan III exhaust plume. 

These correlations show that increasing the eddy viscosity coefficient by a 
factor of 5 over that used in the free-flight case gives a reasonably good 
representation of the measured spreading rate for the 100 K above ambient 
isotherm. 

The principal afterburning species affected by the water injection and 
increased turbulence were found to be CI 2 , CO, and NO. The combined effects 
of increased turbulence and injected water flow twice that of the mass flow 
through the nozzle exit plane are to suppress NO production almost completely, 
to lower the amount of CI 2 produced in the plume over free-flight conditions, 
and correspondingly increase slightly the HCl production. In addition, these 
effects will produce larger amounts of CO than in the free-flight case, but 
the ratio of CO 2 to CO will still be on the order of a hundred to one. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
December 2, 1977 
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TABLE I.- RESULTANT e-FOLDING DROPLETS DECAY TIMES AS A FUNCTION 
OF KINETIC RATE COEFFICIENT 
[^T is assumed to be 373 K]] 


Assumed kf(373 K) 10-20 iq-17 iq-16 iq-IS 

Resultant e-folding 

time, sec 4 x 10“^ 4 x 10“^ 4 x 10“^ 4 x 10“® 4 x lO”"^ 


TABLE II.- DROPLET EVAPORATION TIMES (4 e-FOLDING TIMES) AS 
PREDICTED BY THE DIFFUSION - HEAT FLUX MODEL 


Initial droplet radius, ym 
Evaporation times, sec 


100 

5 


10 1 0.1 

5 X 10-2 5 10-4 5 X 10“6 


J 
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TABLE III.- CHEMICAL SCHEME® 


( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 

(7) 

( 8 ) 

(9) 

( 10 ) 

( 11 ) 

( 12 ) 

(13) 

(14) 

(15) 

(16) 

(17) 

( 18 ) 

(19) 

( 20 ) 
( 21 ) 
( 22 ) 


Reaction 

H + Cl + M = HCl + M 

H + HCl = Cl + H2 

HCl + OH = H2O + Cl 

H+H+MBH 2 +M 

H + OH + M = H 2 O + M 

OH + OH = H 2 O + 0 

OH + H 2 = H 2 O + H 

CO + OH = CO 2 + H 

0 + 0+ M = 02+M 

CO + 0 + M = CO2 + M 

0 + H2 B OH + H 

H + O2 = OH + 0 

0+H+ M=OH+M 

0 + N2 = NO + N 

N + O2 = NO + 0 

NO + 0 + M = NO2 + M 

NO + OH = NO2 + H 

Cl + OH = HCl + 0 

Cl + Cl + M = CI2 + M 

Liquid water + N2 H2O + N2 

N + N+ M = N2+M 

H + CL2 n HCl + Cl 


Rate coefficient 


8.8 

7.2 

1 

1 


!.8 
!.0 
i .6 


X 10- 
X 10- 

X 10' 
X 10- 
10 ' 
10 - 
10- 
10- 
X 10- 
X 10" 
X 10- 
X 10- 
X 10- 


2.3 
1.1 
4.17 X 10 


X 10- 
X 10- 


10 - 

10 - 

10 - 


See text 
5.6 E-30 t 


29x-1 

exp(-4 622/RgT) 

12 exp (-3 250/SjjT) 

29t-1 ® 

28-p-l 

11 exp(-1 000/RgT) 

11 exp (-5 500/RgT) 

13 exp(-600/R-T) 

29t-1 ® 

^^T-1 exp (-4 000/RgT) 

11 exp(-8 200/RgT) 

1® exp (-16 500 /RjjT) 

29t-1 ® 

10 exp (-75 400/RgT) 
I^T exp(-6 250/RgT) 

33 exp(+1 860/RgT) 

12 exp(-30 000/RgT) 

11 exp (-9 900/RgT) 

33 exp(+1 600/RgT) 


-1 


!.0 E 


-10 


®Rate coefficient units are in terms of molecules, cm3, sec. AlCl, AICI2, 
and AI2O3 are included in the plume as inert species. These rate coefficients 
are taken from reference 2. 
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TABLE IV.- SPACE SHUTTLE SOLID- ROCKET MOTOR EXIT 


PLANE AVERAGE VALUES 

[Space shuttle data obtained from Ben Shackleford, 
Marshall Space Flight Center, Mar. 1974] 


Exit plane 
species 

Mole fractions 

AlCl 

0.00003 

AICI 2 

.00008 

AICI 3 

.00001 

AlOCl 

.00002 

AI 2 O 3 

.07975 

CO 

.23205 

CO 2 

.02120 

Cl 

.00225 

Fe 

.0001 1 

FeCl 

.00001 

FeCl2 

.00122 

H 

.00556 

HCl 

.15442 


Exit plane 
species 

Mole fractions 

H2 

0.27812 

H 2 O 

.14068 

NO 

.00001 

N2 

.08401 

0 

.00001 

OH 

.00046 

CIO 

8.1 X 10-9 

CI 2 

6.9 X 10-7 

HO 2 

1.2 X lO-'IO 

NO 2 

5.8 X 10-12 

N 2 O 

8.6 X 10-11 

O 2 

2.0 X 10-7 

N 

0.0 


Exit plane conditions based on an equilibrium model: 


Pexit> 1.00 

'^exit» ^ 2308 

Mexit 2.779 

rexit> 1.810 

Uexit> 2433.8 


TABLE V.- PLUME COMPOSITION AT 1 km DOWNSTREAM OF EXIT PLANE 


[kf(T) = 


Entrained water flow/Original 

gas flow at exit planes ... 2.4 1.3 0.85 0.65 


Species 

Quantity of species, g/m 

HCl 

23 146 

22 206 

21 831 

21 405 

H2O 

271 081 

166 145 

110 054 

85 685 

CO 

360 

97 

76 

69 

CO2 

46 446 

47 200 

47 461 

47 245 

AlpO? 

34 290 

34 719 

34 684 

34 452 

NO 

20 

628 

1 207 

1 805 

CI2 

977 

1 990 

2 457 

2 689 

Liquid water 

63 162 

15 100 

5 072 

4 218 


TABLE VI.- PLUME COMPOSITION AT 1 km DOWNSTREAM OF EXIT PLANE 
AS A FUNCTION OF REACTION RATE COEFFICIENT 
FOR EXCESS WATER FLOW CASE 


jljiquid water + N2 


H2O + N2] 


kf(T) 

|t 


10“20 10“’'7 10 “’'^ 10“''5 


Species 

Quantity of species, g/m 


HCl 

22 215 

23 361 

23 146 

23 701 

CO 

157 

303 

360 

2 307 

CO2 

47 575 

47 238 

46 446 

42 908 

NO 

171 

16 

20 

3 

CI2 

2 144 

1 014 

977 

3 

AI2O3 

34 778 

34 826 

34 290 

33 936 
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TABLE VII.- FINAL PLUME COMPOSITION AS A FUNCTION OF TURBULENT MIXING 


[kf(T) = 10“^^'1 t; this table contains estimates of 
asymptotic values which will be slightly higher 
than the values of 1 km downstream^ 


a 


5 10 20 100 


Species 


Quantity of species, g/m 


HCl 

31 116 

30 608 

30 320 

29 807 

28 649 

Cl2 

1 321 

1 834 

2 030 

2 448 

3 566 

Cl 

27 

37 

43 

65 

174 

CO 

487 

547 

554 

615 

974 

CO2 

62 835 

64 879 

62 670 

62 373 

62 054 

AIoOq 

46 390 

46 390 

46 390 

46 390 

46 390 

NO 

27 

18 

17 

16 

15 

H2O 

366 628 

397 779 

389 870 

392 607 

398 212 

Liquid water . . . • 

85 450 

80 854 

65 163 

60 811 

57 082 


TABLE VIII.- QUANTITIES OF SPECIES FOR FREE FLIGHT CASE 
AND EXCESS WATER CASE WITH a = 5 FOR SAME VELOCITY 
AND NOZZLE EXIT PLANE CONDITIONS 



Free flight 

Combined effects 

HCl, g/m 

29 156 

30 608 

CI2, g/m ^ 

3 282 

1 834 

NO, g/m 

2 029 

18 

CO , g/m 

106 

547 
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Figure 1.- Pad sound suppression scheme. 
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No water in atmosphere 
Water vapor in atmosphere 
Liquid water in atmosphere 


Figure 2.- 


Distance downstream of exit plan, m 

Quenching effects of water on nonreacting shuttle-like plume, 
of water used is very large and identical in both cases. 


Amount 



Center -line temperature 



Distance downstream of nozzle exit plane, m 


Figure 3.- Plume center-line temperature as a function of distance downstream of nozzle 
exit plane for various amounts of water entrained. Evaporation rate coefficient used 
is 10-1°^. 


ru 

uo 


Mole fraction of species along center line 


10 "^ 


10-3 



Distance downstream of exit plane, m 

Figure 4.- Mole fractions of combustion products along plume 
center line as a function of distance downstj^eam of nozzle 
exit plane. Rate coefficient used is Excess flow 

case; Mass HpO/Mass exhaust gas, 2. 


Mole fraction of species along plume center line 



Distance downstream of nozzle exit plane, m 

Figure 5.- Mole fractions of combustion products along plume 
center line as a function of distance downstream of nozzle 
exit plane. Rate coefficient used for water evaporation 
is Mass HpO/Mass exhaust gas, 2. 


Quantity of species in g/m 


Distance downstream of exit plane, m 


Figure 6.- Quantities of chlorine containing species 
(g/m) in plume as a function of distance dovmstream 
of nozzle exit plane. The evaporation rate coefficient 
used is 10“^^\|T. Mass HpO/Mass exhaust products, 2. 




Quantity of species in g/m 



Center -line temperature, 



Distance downstream of nozzle exit plane, m 

Figure 8.- Center-line temperature as a function of distance downstream 
of nozzle exit plane in excess water case for several evaporation 
rate coefficients. Mass H20/Mass exhaust gases, 2 



Quantity of species in g/m 



Figure 9 .- Quantities of chlorine containing species in plume (g/m) 
as a function of distance. Downstream of nozzle exit plane for 
various evaporation rate coefficients. Excess water case 
(Mass H20/Mass exhaust = 2) is considered. 



[)0 800 1000 
e exit plane, m 


reduced by using three 
efficients. 




< 3480 m" > 

L-77-388 

Figure 12.- Infrared scanner photograph of heat pattern from Titan III, 
March 14, 1976. Altitude approximately 3 km. Ambient clouds and 
ground form background intensities. 
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Figure 13.- Infrared scanner photograph of Titan III, May 20, 1975. 
Altitude approximately 1.1 km. 100 K above ambient (400 K) 
isotherm is shown for plume. 
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a = 80 m 

b = 200 m 
c = 280 m 


(a) 400 K isotherm of plume in ordinary coordinates. 

e = 0.5. 


a 



a = 80 m 

b = 200 m 
c = 280 m 


(b) 400 K isotherm of same plume (idealized in 

downstream coordinate system) . 

Figure 14.- 400 K isotherms. 
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Downstream distance, m 


Figure 15.- 400 K isotherm for various turbulence parameters a 
in three cases until downstream distance reaches 80 m. R = 
Mass of entrained HoO/Mass exhaust, 0.25. 




Center -line temperature, 






1. Report No. 

NASA TP-1111 


4 . Title and Subtitle 


2. Government Accession No. 


EFFECTS OF ENTRAINED WATER AND STRONG TURBULENCE ON 
AFTERBURNING WITHIN SOLID ROCKET MOTOR PLUMES 


3. Recipient's Catalog No. 

5. Report Date 

January 1978 

6. Performing Organization Code 


7. Author(s) 


8. Performing Organization Report No. 


Richard I, Gomherg and Richard G. Wilmoth 

9. Performing Organization Name and Address 

NASA Langley Research Center 
Hampton, VA 23665 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, DC 205^6 

15. Supplementary Notes 


L-118T5 


10 . Work Unit No. 


506-21-30-01 

11. Contract or Grant No. 

13. Type of Report and Period Covered 

Technical Paper 

14. Sponsoring Agency Code 


16. Abstract 


During the first few seconds of the space shuttle trajectory, the solid rocket 
boosters will be in the proximity of the launch pad. Because of the launch pad 
structures and the surface of the Earth, the turbulent mixing experienced by the 
exhaust gases will be greatly increased over that for the free-f light situation. 

In addition, a system will be present, designed to protect the lifting vehicle 
from launch structure vibrations, which will inject large quantities of liquid 
water into the hot plume. This report studies the effects of these two phenomena 
on the temperatures, chemical composition, and flow field present in the after- 
burning solid rocket motor exhaust plumes of the space shuttle. The study includes 
results from both a computational model of the afterburning and supporting measure- 
ments from Titan III exhaust plumes taken at Kennedy Space Center with infrared 
scanning radiometers. 


17. Key Words (Suggested by Author (s)) 


18. Distribution Statement 


Afterburning 

Space shuttle 

Launch vehicle effluents 


Unclassified - Unlimited 


Subject Category 


19. Security Oassif. (of this report) 1 

20. Security Classif. (of this page) , 

1 21. No. of Pages | 

1 22. Price* 

Unclassified j 

Unclassified 

1 

$^.50 


* For sale by the National Technical Information Service, Springfield. Virginia 22161 


IVASA-LangTey, 1978 



THIRD-CLASS BULK RATE 


National Aeronautics and 
Space Administration 

Washington, D.C. 

20546 

Official Business 

Penalty for Private Use, $300 


Postage and Fees Paid 
National Aeronautics and 
Space Administration 
NASA-451 



4 1 1U,A, 121277 S00903DS 

DEPT OF THE AIR FORCE 
AF WEAPONS LABORATORY 
ATTN: TECHNICAL LIBRARY (SOL) 
KIRTLAND AFB NH 87117 




POSTMASTER: 


If Undeliverable (Section 158 
Postal Manual) Do Not Return 


V 




